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Control of the Josephson current by varying a gate current has recently been demonstrated in both 
4-terminal and 3-terminal junctions. We show that, when the the gates are weakly coupled to the 
Josephson junction, the Josephson current versus gate current (or versus gate voltage) relation is 
the same for both the 4- and 3- terminal geometries. At low temperature, the supercurrent switches 
abruptly as a function of the gate voltage, but only slowly as a function of the gate current. 
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I. INTRODUCTION 



Transistors with superconducting sources and drains 
have been fabricated over the last 10|years, using field ef- 
fect control of the supercurrent flow.El For the field effect 
to change carrier density typically requires a gate voltage 
of the order of volts, while superconducting energy gaps 
{Egap — 2A) are at best of order millivolts. To over- 
come this incompatibility in voltage scales, van Houten 
suggested using quantum confinementB in the channel to 
incorporate a new small energy scale into the transistor. 
Due to the added confinement energy, a small change in 
gate voltage can produce a large change in the supercur- 
rent flow. 

A small energy scale occurs naturally in the supercon- 
ducting transistor, even without the added quantum con- 
finement. The two superconductors themselves form a 
special type of 'barrier' for electrons, where the electrons 
'Andreev reflect' for the superconductor as a hole.l3'cl Re- 
peated Andreev reflections from the two superconductors 
form a closed path, leading to energy level quantization 
inside the device channel. These 'Andreev energy levels' 
form inside the superconducting gap, with more bound 
levels forming as the channel length increases (just as in 
the standard type of quantum well). A large fraction 
of ihe Josephson current flows through the Andreev lev- 
elsQ 13, where a quasi-particle moving once around this 
closed path corresponds to the transfer of a pair of elec- 
trons across the device channel. 

Figure |l| shows a type of gated Joseohson junctions 
used in two different recent experimentsBB to control the 
Josephson effect. The gate terminal are all normal metals 
(Ni,N2) held at voltages (Vi,V2)- The superconductors 
(S) are grounded. Electrons can tunnel from the gates 
into the normal region (N) of the Josephson junction with 
probabilities {Ti,T2). Setting the tunneling probability 
T2 = decouples one of the gates from the Josephson 
junction, forming a^three terminal junction. 

Morpurgo et aLEI have recently used current flow 
through the two opposing gates to switch off the Joseph- 



son current in a 4-terminal junction. Control of the 
Josephson current with a gate voltage (or current) in 
a more convcntional_3-terminal geometry has also been 
shown by SchalpersLl et al.. Both recent experiments 
are based on the mechanism of van Wees et al.Q that, 
because the superconductor cannot inject electrons into 
the bound Andreev levels, the electrochemical potential 
/is of the bound levels will float to the probe voltage 
(e.g. hb — eVi in a 3-terminal geometry) rather than 



to the electrochemical 



potential of the superconductors 
have also used this principle to 
analyze the the current flow in different types of gated 
Josephson junctions. 
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FIG. 1. A four terminal Josephson junction. The gate 
currents /qi and /g2 passing through the normal region con- 
trol the supercurrent I by altering the occupation of bound 
energy levels in the normal region. 



II. FOUR TERMINAL JUNCTION 



Consider next the Andreev level occupation in a 4- 
terminal Josephson junction. The occupation factor fs 
for the Andreev levels (inside the normal region (N) of 



1 



the Josephson junction where the Andreev levels exist) 
is now given by 



T1+T2 



(1) 



The occupation factors /i and /2 in Eq. arc the Fermi 
factors of the two gates, namely /i — f{E — eVi) and 
/i = f{E — eV2). Here Vi and V2 are the gate voltages 
and f{E) = l/[l + exp{E/kBT)] is the equilibrium Fermi 
factor. Equation ([l]) states that the occupation factor for 
an electron in the normal region is simply the probability 
that the electrode originated from probe i — \,2 times 
the occupation factor fi of probe i. The occupation fac- 
tor in Eq. (|^) is the same as for a small metallic grain, 
semiconductor quantum dot, etc. connected to two nor- 
mal metal leads. For the 3-terminal device {T2 — 0), 
Eq. gives the occupation factor of the bound levels is 
simply the Fermi factor of the gate /s = /i . The bound 
level occupation factors in the 3- and 4- terminal devices 
are shown in Fig. ^. 
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FIG. 2. Occupation factor /s(-B) of the supercurrent car- 
rying energy levels {\E\ < A) versus electron energy (refer- 
enced to the Fermi level of the superconductors fis = 0) in the 
(a) three terminal and (b) four terminal Josephson junctions. 
Both occupation factors give the same terminal currents, even 
though the occupation factor in (b) is clearly not an equilib- 
rium Fermi factor. Heating the distributions (dotted lines) 
smears the switching of the Josephson current, so that the 
device performs better with cold electrons (solid). 



the same for both 3-terminal and 4-terminal devices. The 
'bound level' current Ig flows in the energy range-within 
|A| of the Fermi level {\E\ < A), and is given bytll 



Ib — ^ In{En)fB{En) 



(2) 



Here En is the energy of the Andreev level and /„ is the 
current carried by the level. In the weak coupling limit 
(Ti ^ 1, T2 ^ 1), where the wavefunctions inside the 
normal region are essentially unchanged by adding the 
gates, /„ is the current carried by bound levels in the 
Josephson junction before adding the gates. Inserting 
Eq. d) into Eq. (|) gives 



n 

rr _,rr. Y.^n{En)f{En-eV2). (3) 
1 + ^2 — 



TI + T2 
T2 



Equations (|l|)-(|^) are easily generalized to any number of 
gates connected to the Josephson junction. 

Morpurgo et al. used a symmetrical device (for which 
we define the gate coupling Ti =T2 = e), passing a gate 
current Iqi = —Ig2, so that the, gate voltages satisfy 
Vi = —¥2- One can easily verifytJ that the bound level 
current Ib does not depend on the sign of the gate voltage 
in a 3-terminal structure, so in this limit Eq. (^) reduces 
to 



/-(l-gato) 



(2 — gates) 
B 



(4) 



To plot the Josephson current versus gate current (for 
the symmetrical 4-terminal device in the weak coupling 
limit), we can ther efore use the same formalism previ- 
ously develop edlilE3 for the 3-terminal junction. Specif- 
ically, we use the one-dimensional Josephson junction 
model of Ref. |l^, where the normal region extends over 
Q < X < L and includes a tunnel barrier having trans- 
mission probability T at a; = a. Ref. ^ uses the sim- 
plest type of scattering matrix to describe the connec- 
tion of the gate to the Josephson junction, which neglects 
backscattering of quasi-particles due to the gate. Addi- 
tional physics may develop if one allows the gate to al 
the quasi-particle trajectories through backscattering. 



The Josephson current flow through the either desjice 
in Fig. is divided into two different energy regional as 
I ^ Ib + Ic- The 'continuum' current Ic flows at ener- 
gies outside the superconducting gap {\E\ > A), and is 
essentially not aiiected by the additional gates (neglect- 
ing gate leakage) .llll The current Ic is therefore essentially 
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FIG. 3. Terminal characteristics for the short (L <C ^o) Josephson junction, (a) Differential conductance along the gate 
versus gate voltage, (b) Josephson current versus gate voltage, (c) Gate current versus gate voltage, (d) Josephson current 
versus gate current. The paramater T is the probability to transmit across the normal region. 



Even though the 4-terniinal device has a 'double step' 
Fermi distribution, shown in Fig. ^(b), this unusual 
nonequilibriuni distribution function is irrelevent to the 
Josephson current switching. The current is the same as 
if the bound levels had the standard equilibrium Fermi 
distribution of Fig. ||(a). The temperature of the elec- 
trons in this distribution, whether they are 'hot' or 'cold', 
is also not important to the switching effect. Switching 
of the Josephson current versus gate voltage or gate cur- 
rent will be sharper if the electrons remain 'cold'. For 
the 3-terminal device having the distribution function of 
Fig. ^(a), there is also a voltage drop of Vq across the 
NS junctions. The total voltage drop between the two su- 
perconductors is still zero, as the voltage difference across 
the left NS junction is the negative of the voltage drop 
across the right NS junction. In the 4-terminal device, 
where the average electrochemical potential is equal to 
that of the superconductors (^s = 0), there is no voltage 
drop across either NS junction. The presence of absence 
of these voltage drops does not (to lowest order) affect 



the supercurrcnt flow. 

A. Short Junction 



Figure ^ shows the terminal characteristics of a short 
(L <C ^Oj where L is the length of the normal region and 
^0 the coherence length of the superconductors) Joseph- 
son junction. The gate couplings in Fig. ^ are weak, with 
T\ — Ti — \% = We fix the superconducting phase 
difference at </) = .Gtt. The differential conductance along 
the gate shows a resonance whenever the gate voltage 
crosses a new Andreev level in Fig. ^(a), a tool which 
can be used to perform 'Andreev level spectroscopy'. The 
gate current versus voltage in Fig. ||(c) is simply the inte- 
gral of Fig. ||(a), showing steps in the gate current instead 
of peaks in differential conductance. The Josephson cur- 
rent versus gate voltage in Fig. ^(b) switches abruptly to 



3 



zero whenever a new Andreev level is populated. Reduc- 
ing transmission T across the normal region from 100% 
to 12% in Fig. ^(b) suppresses the Josephson current and 
pushes the Andreev level closer to the gap edge, pushing 
the switching voltage close to eVc — A. 

The Josephson current versus gate current is shown 
in Fig. ^(d). Both the Josephson current switching in 
Fig. |3|(b) and the increase in the gate current in Fig. ^(c) 
occur over over the same range of voltage, namely the 
width (in energy) of the Andreev level. The Josephson 
current therefore decreases quite slowly and gradually 
with increasing gate current. The fractional occupation 
of the single Andreev level in this junction {L <C ^o) also 
increases slowly from /s = to = 1 during the switch- 
ing event. Contrast the sharp switching characteristics in 
the Josephson current versus gate voltage from Fig. ||(b) 
(roughly a square wave) to the relatively slow change of 
Josephson current versus gate current in Fig. ^(d) (ap- 
proximately triangular). 
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FIG. 4. The differential conductance along the gate versus 
gate voltage shows more Andreev levels for the long [L ^ ^o) 
Josephson junction. The long junction is either (a) ballis- 
tic, (b) an 'asymmetric' tunnel junction, or (c) a 'symmetric' 
tunnel junction. 
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FIG. 5. For a long Josephson junction, which is either bal- 
listic (solid), a 'symmetric' tunnel junction (dashed), or an 
'asymmetric' tunnel junction (dot-dashed), we show (a) the 
Josephson current versus gate voltage, (b) the gate current 
versus gate voltage, and (c) the Josephson current versus gate 
current. 



B. Long Junction 



Figures ^ and || show the same terminal characteris- 
tics for a long {L 3> Co) Josephson junction. In long 
SNS junctions the position of the tunnel barrier has 



some influence on the size of the Josephson currentQ'talla. 
We therefore consider a 'ballistic' junction with T — 1, 
an 'asymmetric' junction with the tunnel barrier 1/5 of 
the distance across the normal region (a = L/5 and 
T = 12%), and a 'symmetric' junction with the tunnel 
barrier in the middle of the normal region {a = L/2 and 
T = 12%). The gate is again only weakly coupled to the 
Josephson junction, with Ti = T2 = 1% = e. We again 
fix = .67r. This particular long junction has three An- 
dreev levels in the energy range within A of the Fermi 
level (since L = 6.6^0), as can be seen from the differ- 
ential conductance along the gate versus gate voltage in 
Fig. ^. Only small differences appear in the differential 
conductance along the gate between the different types 
of long junctions in Fig. ^. 

Figure ^ shows (a) the Josephson current versus gate 
voltage, (b) the gate current versus gate voltage (the 
integral of Fig. and (c) the Josephson current ver- 
sus gate current for the three different types of long 
Josephson junctions. The Josephson current in Figs. ||(a) 
and (c) changes sign whenever the gate voltage is such 
that pi|i*ipw Andreev level becomes populatodjps depopu- 
latedEHlil (sometimes called a '7r-phase shift 'll3llJ). A sign 
change of the Josephson current only occurs in long junc- 
tions, as filling or emptying the single Andreev level in a 
short junction forces the Josephson current to zero. The 
Josephson current remaining in Fig. |^(a) and (c) when 
all the Andreev levels are filled {eVa > A), is the 'con- 
tinuum' current Ic- 

The Josephson current switches abruptly with gate 
voltage in Fig. ||(a), but only gradually with gate cur- 
rent in Fig. ^c). The with in voltage over which these 
currents change is again the Andreev energy level width. 
Andreev level is filling changes the gate current and 
Josephson current simultaneously, leading again to only 
slow variation of Josephson current with gate current in 
Fig. |^(c). The same reasoning explains why the shape 
of Fig. ||(a) is roughly a square wave, while Fig. ||(c) is 
roughly triangular. 

The ballistic junction (solid lines) always has the 
largest Josephson current, as seen in Fig. H(a) and (c). 
The symmetric junction (dashedj-linfis) displays the so- 
called 'giant' Josephson currenlll3'E3 as the switching 
event near eVc — .35 A in Fig. ^(a), but the asymmetric 
junction can also display switching currents of the same 
magnitude. One indeed expects the symmetric tunnel 
junction to have laxge switching currents compared to its 
Josephson currenttJ, but the asymmetric junctionj-can 
have these as well when the Andreev gaps are smallEj. 

Since Ref. || uses metals to fabricate the Joseph- 
son junction, the quasi-particle trajectories forming the 
bound Andreev levels will not be .the simple types of scat- 
tering trajecties assumed in RefsOO, but will be altered 
by impurity scattering inside the Josephson junctionJlj'Efl 
In such dirty junctions at very low temperature another 
energy scale becomes important (the Thouless energy, set 
by the time for an electron to diffuse across the Joseph- 
son junction). This extremely low energy scale does not 
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increase as the superconductor improves, and is not es- 
sential to understanding the Josephson current switching 
(or 'TT-phase shifting'). The Josephson current changes 
sign in long junctions whenever a new Andreev level is 
populated or depopulated. Populating a new Andreev 
level in the short junction simply forces the Josephson 
current to switch off. 



C. Maximizing Current Gain 



We now consider how to maximize the change of 
Josephson current over the change in gate current while 
the Josephson current switches. For simplicity we con- 
sider only the ballistic junction, T — 1. Following Ref. |l^, 
the change in Josephson current (A/) is 



(A/) 



(5) 



with £,{En) the energy dependent coherence length. The 
change in gate current near a resonance is 



(A/g) 



(6) 



where Ti — T2 = e is the gate coupling. The ratio of 
Josephson to control current is therefore 



(A/) 
(A/g) 



1 



(7) 



If we could take the gate coupling to zero, the current 
gain would become infinite. However, inelastic or phase 
breaking scattering sets a lower limit on the gate cou- 
pling needed for the gate control to affect the Josephsaa. 
current If the inelastic scattering time is r^, we havetll 



VFT4, 

The maximum current gain then becomes 



(A/) 



(A/g) 



(8) 



(9) 



where the arrow indicates the short junction limit. The 
maximum current gain in the ballistic short junction limit 
is therefore set by the ratio of the superconducting energy 
gap to the level broadening induced by inelastic or phase 
breaking scattering. The device of Ref. in the short, 
ballistic junction limit, shows a current gain of about 20 
for small gate currents. 



III. CONCLUSIONS 



We have shown the dependence of the supercurrent 
flow in a three or four terminal Josephson junction on the 
current passing through ajMrmal metal gate, motivated 
by two recent experimentsMEl When the gates or gates are 
weakly coupled to the Josephson junction, desirable for 
better isolation of gate (control) current from the Joseph- 
son current, there is little (or no) difference between the 
terminal characteristics of the three versus four terminal 
geometries. The Josephson current in either geometry is 
a slowly varying function of the gate current, even for 
the idealized one-dimensional model we have considered 
in this paper. The Josephson current varies much more 
rapidly as a function of the gate voltage. Measuring the 
differential conductance along the gate would also permit 
direct observation (spectroscopy) of the current carrying 
energy levels inside the superconducting gap. _ 

In both the transistor farbricated using metalsEl and in 
compound semiconductorga, the gates are not well iso- 
lated from the device channel. Instead of the tunnel bar- 
riers drawn in Fig. |l], only impurity scattering inside the 
Au metalB or two-dimensional electron gasO, along with 
possible contact resistances, provide partial gate isola- 
tion. Since quasi-particles which are normally trapped 
inside the Josephson junction can now escape by leaking 
into the gate, the bound Andreev levels in either device 
will therefore be very broad. Broad Andreev levels not 
only broaden the Josephson current switching, but pro- 
duces gate currents in experiments which are compara- 
ble to the supercurrents. As pointed out in Ref. ||, The 
two dimensional (planar) geometry of the devices used in 
the two exjieriments also introduce a type of geometrical 
broadenindlj which needs to be studied in this context. 
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